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ABSTRACT 

Measurements  of  sea  clutter  using  liigh- resolution  radar  indicate 
that  the  cluiter-eross-seetion  returns  lollow  a  log-normal  probability 
density  function  more  closely  than  the  usually  assumed  Rayleigh  law. 
This  report  develops  the  theory  for  the  detection  of  a  steady  signal  In 
log-normal  clutter  by  first  using  a  single  pulse  and  then  by  using  the 
cum  of  /V  pulses  integrated  noneoherently .  Plots  of  the  probability  den¬ 
sity  of  the  envelope1  of  the  signal  plus  clutter  show  the  function  to  be 
blmodnl,  an  unexpected  result.  Curves  are  presented  for  the  threshold 
bias,  normalized  to  the  median  clutter  voltage,  versus  the  probability 
of  false  alarm  for  several  values  of  the  standard  deviation  ,i  and  for 
various  values  of  /v.  Probability  of  detection  curves  arc  preseated  for 
«  -  3,  0,  and  9  dll,  for  1,  3,  10,  and  30  pulses,  and  for  false  alarm 
probabilities  from  10  1  to  10  H.  The  ratio  of  signal  to  median  clutter 
required  for  detection  increases  markedly  as  ,>  increases  because  of 
the  highly  skewed  clutter  density. 
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SUMMARY 

Sen  rhiHee  mcnmiremeiUH  made  by  NHl,  union  a  high  icMtluli<m  radar  show  that  I lw* 
chiller  cross  him  (tun  tends  to  be  big  normally  distributed,  A  idudy  o|  tin  by*  mutual 
clutter  model,  theretoiv,  seemed  appi  update.  ami  this  report  present*  tlu  ro*ult*  »>| 
such  a  study, 

V'li’Hl,  a  derivation  in  presented  o|  i|»>  Ion- normal  density  lumtlon  required  to  e*tnh 
lisl\  a  family  of  threshold  values  (or  lalse  alarm  proltuhllitles  varying  from  10  to  10 
uml  for  selected  values  of  standard  deviation  ..  Riots  of  the  flutter  density  show  highly 
shewed  tall,-,  as  .•  tnere.iset  result  out  la  ureatly  Increased  Idas  values  as  eomiiarnl  to  the 
Rayleigh  motlel, 

Next,  the  distribution  of  the  envelope  ol  signal  plus  flutter  lor  a  simile  pulse  is  de 
rived,  Plots  of  probability  densities  show  a  Id  modal  function,  which  Is  an  unexpected 
result,  The  btmodal  shape  holds  tor  all  values  of  signal- to  noise  ratio  tS  N)  ami  for  all 
values  ot  but  the  dual  peaks  draw  closer  together  and  the  dtp  In  iomcs  less  pronounced 
as  ■  incre.iN  >s,  Formulas  lor  the  prohabilttv  id  deteethm  are  given,  hut  these  could  not 
he  evaluated  In  closed  hu  m,  so  they  were  programmed  lor  the  NHl.  CDC  3M00  computer, 

Then  the  probabilities  ol  false  alarm  and  detection  are  developed  for  tire  sum  ot  \ 
pulses  integrated  noncoherent ly.  The  general  method  used  is  thul  of  the  characteristic 
function,  which  essentially  requires  two  .ipplleatious  ol  tl>c  Fourier  integral,  Here  again 
no  solution  in  closed  form  seemed  readily  obtainable,  so  tire  t'l)t'  31UH)  was  imrgrammed 
for  numerical  evaluation,  using  the  fast  Fourier  transform  technique. 

Sets  of  curves  are  presented  for  the  probability  of  detection  fur  values  of  ..  3,  tl,  and 

9  dll,  for  A  1,  3,  10,  and  30  pulses,  and  for  false  alarm  prolmbt  titles  from  10  ’  to  10  h. 
The  y  N  required  for  defection  for  ■  3  dll  are  about  the  same  as  those  for  the  Rayleigh 

model,  as  would  bo  expected  since  the  log-normal  tail  for  tilts  ease  compares  to  the 
Rayleigh  tall.  However,  as  ••  increases,  the  required  S  N  also  rises  markedly  as  a  result 
of  the  Increasingly  skewed  distributions.  It  should  be  noted  Hull  the  y  N  defined  In  lids 
report  uses  the  median  clutter  value  for  noise  and  not  the  usual  vms  value. 


INTRODUCTION 

Clutter -cross -sect Ion  measurements  of  sea  return  taken  with  high- resolution  radar 
show  more  of  a  tendency  to  follow  a  log-normal  probability  distribution  than  the  Rayleigh 
distribution  usually  assumed.  Data  reported  by  NRL  (1'  using  an  X-band  radar  with  a 
0, 02-n  see  pulse  are  plotted  in  Fig.  1  on  probability  paper  with  the  clutter  cross  section 
<»,  in  dB  on  an  arbitrary  scale.  A  log-normal  distribution  Is  Indicated  by  these  points 
falling  in  a  straight  line.  It  Is  seen  that  the  cumulative  probability  closely  fits  a  line 
between  5 ''('  and  80‘X.  Accurate  data  at  the  extremes  are  difficult  to  obtain,  because  ther¬ 
mal  noise  limits  the  low  ..t  values  and  the  system  dynamic  range  limits  the  high  „r  values. 

Studies  made  by  Miner  (2)  and  Ballard  (3)  also  consider  the  log-normal  description  of 
the  sea-clutter  cross  section  and  relate  the  standard  deviation  of  the  distribution  to  the 


1 


S  .  I  i  it  I  'Hut , 


i  t'M|  M  .  NA%\  '  \  >  I  ii'S  ,v  tj|tt 

F'g,  t  •  t  Vehubilit  y  dlst  r  Unit  sun  lor  lin¬ 
ing  -nm'iiuil  ilutter  u  mile  I ,  w  i  l  h  e  \  pe  r  i  - 
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railin'  llluminulod  patch  area.  Ballard  actually  sot  up  a  lop- normal  probability  model  but 
did  not  calculate  detection  curves.  The  present  report  follows  much  of  the  same  theoreti¬ 
cal  method  and  extends  the  mode)  to  include  the  integration  of  N  pulses. 


UKi-NQHMAl.  CLUTTER  MODEL 

If  we  assume  that  clutter  cross  section  Is  log-  normally  distributed,  then  the  proba¬ 
bility  density  is  given  by 


where  ",  is  the  clutter  cross  section,  ,.;II  is  the  median  value  of " and  ,,  is  the  standard 
deviation  of  in,,  .(natural  logarithm).  Since  cross-section  values  are  measures  of  power, 
the  standard  deviation  (s.d.)  is  expressed  in  "natural  units”  (4)  by  the  relationship 

"(natural  units)  (0.1  |n  10)  ,.(,115)  0. 2303  ,.(<IB) .  (2) 

'Hte  value  of  in  Fig.  1  is  in  dB,  and  „  must  be  converted  into  natural  units  by  Eq.  (2)  for 
use  in  Eq.  (1). 

In  Fig.  l  the  log-normal  probability  distribution  function  is  plotted  for  ,,  -  3  dB  and 
,,  6  dB.  The  dashed  line  is  a  close  fit  to  the  NRL  experimental  results  on  sea  clutter 

previously  mentioned.  For  this  line.  »  5.4  dB.  From  Ref.  3,  a  value  of  =  4.7  dB 

resulted  from  a  clutter  patch  size  of  3840  sq  yd  in  a  sea  state  *t.  From  the  NRL  studies 
(1),  ,  ■■  5.4  dB  corresponds  roughly  to  a  clutter  patch  size  of  135  sq  yd  in  sea  state  1.* 

To  determine  the  probability  of  false  alarm  and  subsequently  the  probability  density 
function  of  a  signal  plus  clutter,  Eq.  (1)  must  first  be  transformed  into  a  voltage-amplitude 


In  general  >  increases  directly  with  sea  state  and  inversely  with  clutter  patch  size. 
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density.  II  we  lot  v  bo  the  envelope  of  the  radar  echo  from  the  clutter  area  ,  then 
\  .  1  *’  or  „r  hv-  .  Using  this  substitution  In  Eq.  (1)  yields 


where  vm  is  the  median  value  of  v,.  and  .■  remains  the  s.d.  of  in<>(, ,  Then,  .<  is  still  mea¬ 
sured  according  to  Eq.  (2). 


The  probability  of  false  alarm  for  a  voltage  bias  level  /( is  given  by 


r-  i 


c\p  ( - r 2)  dt . 


where  ;  /(  v.i(  is  the  normalized  bias  level.  The  integral  in  Eq.  (4)  is  related  to  the  error 

function,  which  is  well  tabulated. 


DETECTION  OF  A  NON  FLUCTUATING  SIGNAL  IN  CLUTTER 

Let  us  assume  a  radar  pulse  return  from  a  target  which  is  a  steady  signal  of  constant 
amplitude  v  .  The  clutter  return  will  be  assumed  to  have  the  log-normal  amplitude  density 
function  as  given  in  Eq.  (3)  and  a  phase  that  is  uniformly  random  over  -.7  .  The 

analysis  can  be  performed  as  though  the  radar  were  operating  as  a  cw  system.  If  the 
signal  and  clutter  voltages  add  linearly  prior  to  envelope  detection,  then  we  can  assume 
that  the  signal  plus  clutter  may  be  represented  by  the  two  components 

X  vs  •  ms  </j(.  (5a) 

and 

1'  vk.  sm  <.'v.  (5b) 

where  x  and  Y  are  the  in-phase  and  quadrature  components  of  the  signal  plus  clutter. 

We  can  define  two  new  variables  x  and  Y'  such  that  their  joint  probability  is  that  of 
the  clutter  alone: 

X'  X  -  Vs  V,  cos  </v  (6a) 

and 

Y‘  Y  r  sill  ,.'v.  (6b) 

Then,  the  relationship  between  these  joint  distributions  is 

>w  <>r 

itv  .  0  </i 

pi  v  ,<!>r )  ilv  .di!>  1  1 

f>. '  r nr 

ih/i 


[AX'  ,Y' )  cU'tlY'  v(  p(X'.K')  clX'clY' . 


(7) 
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But  since  vc  and  <t>c  are  independent  variables,  p(  vc,<t>c)dvcd<j>c  .  P(  vt.)(/vrp(./v  )d</v 
-  p\.vc)dvc{l/2n)ckf>c,  lienee,  from  Eqs.  (3)  and  (7)  we  have 


_  JL  exo 

- 1 

■ 

•c 

;  r. 

L  2(7  2  \ 

«7vcfA/>t.  ■ 


From  Eqs,  (6)  this  becomes 


P<X,F)  = 


IT  o  UX  -  Vs)2  +  v2l 

Finally,  using  Eqs.  (5)  in  Eq.  (9)  yields 


exp 


2a J 


1  /2  In  ^ 


(8) 


(9) 


p(v,<? i) 


exp 


2 17  a  (v^  -  2v  v  cos  <£  +  v?) 


2  In 


y  v  - 2 vvs  cos  0  » 

V 

V 

m 

/ 

(10) 


where  v2  =  x2  +  v2,  x  -  v  cos  <t>,  Y  *  v  sin  4>,  and  the  transformation  from  Eq,  (9)  to  Eq.  (10) 
is  analogous  to  the  one  used  in  Eq.  (7)  with  0  uniformly  random  from  -n  to  +ff.  Now  v 
represents  the  envelope  amplitude  of  signal  plus  clutter.  The  probability  density  function 
of  v  alone  is  obtained  from  Eq.  (10)  by 

p(v)  2  I  p(v,</>)  d</>.  (11) 

Jo 

Attempts  to  evaluate  Eq.  (11)  in  closed  form  were  not  successful,  so  computer  numerical 
evaluation  was  employed.  However,  before  illustrating  the  results,  it  is  more  meaningful 
to  transform  Eq.  (10)  into  an  expression  containing  a  “psuedor’  S/N,  actually  a  signal-to- 
clutter  ratio,  since  vs  is  not  easily  obtained.  If  we  let  vH/vm  in  Eq.  (10)  the  probability 
density  becomes 


p(x,0 ) 


_ _ X _ _ 

p  V 2 7T  n  (x2  -  2xr  cos  </>  +  r2) 


exp 


--  (in  yfx2  -  2xr  cos  <t>  t  r2)2 
<7  2 


(12) 


where  x  =  v/vm  is  now  the  normalized  envelope  amplitude  of  signal  plus  clutter.  Then  the 
probability  density  function  of  x  is 


p(x)  2  1  p(x,vi)  d<l>.  (13) 

J  0 

The  probability  of  detection  for  a  single  pulse  of  a  steady  signal  in  log-normal  clutter 
for  a  normalized  bias  level  y  is  thus 


r™  r”  ry  rn 

pn  '  I  P(*>  dx  -  2  p(x,<A)  1  -  2  l  p(x ,rf>)  cforfx.  (14) 

In  Fig.  2  plots  of  Eq.  (3)  for  various  values  of  <>  are  compared  with  the  Rayleigh 
distribution  for  which  the  median  value  has  been  set  equal  to  vm  .  It  should  be  noted  that 
the  tail  of  the  Rayleigh  distribution  falls  off  about  as  rapidly  as  that  of  the  log-normal 
for  «=3  dB,  but  the  tail  falls  off  less  rapidly  as  a  becomes  larger  in  the  log-normal 
distribution.  For  these  larger  «  values,  then,  the  bias  values  for  a  fixed  p  are  expected 
to  be  greater  for  the  log-normal  than  for  the  Rayleigh  distribution.  Plots  of  bias  values, 
obtained  from  Eq.  (4),  illustrate  In  Fig.  3  how  rapidly  the  bias  does  increase  as  „  increases. 
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VALUE  OF  THE  CLUTTER  CROSS  SECTION,  vc/*m 


Fig.  2  -  Log-normal  probability  density  for  a  =  3,  6,  and 
9  dB  and  the  Rayleigh  density  tor  the  median  value  v, 


Fig.  3  -  Bias  levels  as  a  function  of  false 
alarm  probability  with  -  3,  6,  9.  and  12 
dB  for  N  -  1,  a  single  pulse  of  log-normal 
clutter 
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To  plo*  the  probability  density  function  in  Eq.  (13)  the  integral  had  tube  evaluated 
numerically,  since  no  closed  form  solution  was  found.  The  CDC  3800  was  programmed 
using  Simpson’s  rule  for  the  integration.  Figure  4  shows  a  family  of  curves  for  the 
probability  density  of  the  envelope  of  signal  plus  clutter  for  »  =  3  dB  and  for  several 
values  of  r  vs  vm  =  S/N  which  has  been  designated  as  a  voltage  S/N.  The  double  hump 
or  bimodal  curve  is  most  unusual  and  surely  not  predicted.  Note  that  this  curve  is  t he 
envelope  distribution  and  hence  assumes  linear  detection.  The  values  of  S/N  shown  arc 
voltage  ratios.  Figure  5  shows  curves  for  =  6  dB.  Again  the  double  hump  is  present 
but  less  pronounced.  In  Fig.  6  by  increasing  »  for  a  fixed  value  of  r  (i.e.,  S/N)  the  dis¬ 
tributions  remain  bimodal  but  the  humps  tend  to  converge  and  the  dip  becomes  less 
pronounced.  Figures  4  and  5  show  that  as  r  increases,  the  curve  shapes  do  not  change 
substantially. 

A  word  of  caution  is  in  order  regarding  the  definition  of  the  signal-to- noise  ratio 
r  =  S/N  =  vs  vm  .  It  should  be  noted  that  this  is  not  the  usual  definition  of  the  signal-to- 
noise  ratio,  where  the  noise  is  taken  to  be  the  rms  value.  Herein  the  noise  (which  is 
clutter)  is  taken  to  be  the  median  value.  This  is  a  convenient  definition  to  use  and,  also, 
one  which  would  have  considerable  engineering  significance,  since  clutter  is  frequently 
expressed  in  terms  of  its  median  value.  However,  in  subsequent  parts  of  this  study 
where  n  pulses  are  integrated,  care  must  be  exercised  in  interpreting  r.  Since  the  log¬ 
normal  model  yields  a  two-parameter  distribution,  its  behavior  is  more  difficult  to 
analyze  and  interpret  than  the  Rayleigh  model. 


DETECTION  USING  THE  SUM  OF  /V  PULSES 

The  straightforward  method  of  determining  the  probability  density  for  the  sum  of  N 
envelope  detected  pulses,  assuming  they  are  independent,  is  by  use  of  the  characteristic 


Fig.  4  -  Probability  densities  of  a  single- pulse  signal 
plus  Log-normal  clutter  for  <>  -  i  dB  and  voltage  S/N 
0,  l,  1,  and  6.  Note  the  double  peaks. 


Fig.  5  -  Probability  densities  of  a  single-pulse  signal 
plus  log-normal  clutter  for  »  =  6dB  and  voltage  S/N 
0,  Z,  -1,  and  6 


normal  ized  envelope  amplitude  of  r.ioNAi  plus  clutter,  »•  e>« 


Fig.  (>  -  Probability  densities  of  a  single-pulse  signal  plus 
log-normal  clutter  for  voltage  S/N  ^  and  ■>  i,  6,  V,  and 
Id  dB. 
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function.  Starting  with  the  probability  density  of  the  envelope  of  signal  plus  clutter  given 
by  Eq.  (13),  the  characteristic  function  is  given  by 

<i>( w)  =  |  p(x)  exp(  twx)  dx,  (15) 

JQ 

which  is  the  Fourier  transform  (FT)  of  p(x),  since  x  0.  The  characteristic  function  for 
the  sum  of  N  independent  pulses  is 


(16) 


It  is  noted  that  in  general  ./>(w)  is  a  complex  function.  Using  the  inverse  FT  of  Eq.  (16) 
gives  the  probability  density  for  the  sum  of  N  pulses: 

P*(x)  ~~  j"  exp(-jwx)  dw.  (17) 

Plots  of  Eq.  (17)  ‘or  N  =  1,  2,  and  >  and  for  S/N  -  2  are  shown  in  Fig.  7.  Note  that  the 
double  hump  disappears  as  H  increases.  Also  the  distribution  spreads  out  as  h  increases 
and  begins  to  lock  more  like  a  normal  distribution. 

For  the  sum  of  pulses,  the  probability  of  false  alarm  is 

P/ytF/l)  j  pN  (x)  c/x  I  -  I  PN(x)  dx,  (18) 

J  y  j  0 

where  r  =  0;  and  the  probability  of  detection  is 

PtfilVt.  )  1  |  p/y(  x)  dx ,  (19) 

J  0 

where  py(x)  is  evaluated  for  the  parameter  r  /  0.  All  attempts  to  find  a  closed-form 
expression  for  the  integral  in  Eqs.  (18)  and  (19)  proved  fruitless.  Fortunately,  the  new 
digital  technique  known  as  the  fast  Fourier  transform  (FFT)  made  it  possible  to  evaluate 
these  integrals  numerically  in  a  reasonable  time. 


PROBABILITY  OF  DETECTION  RESULTS 

Bias  values  for  n  =  3,  10,  and  30  were  determined  from  Eq.  (18)  by  numerical  inte¬ 
gration  for  lalse  alarm  probabilities  from  10 1  to  10"  and  for  «  -  3,  6,  and  9  dB.  These 
are  plotted  in  Figs.  8  through  10.  Curves  for  bias  values  are  repeated  in  Figs.  11  through 
13  with  n  varied  for  a  fixed  «  value. 

Using  these  bias  values,  Eq.  (19)  is  evaluated,  giving  the  probability  of  detection 
versus  S/N  per  pulse,  in  values  of  Vn  from  0.1  to  0.9999,  Figures  14  through  17  are  for 

3  dB.  A  summary  of  PD  for  PFA  -■  10'6  and  N  I,  3,  10,  and  30  is  presented  in  Fig.  18. 
In  Figs.  19  through  23  detection  curves  are  given  for  .>  -  6  dB,  and  in  Figs.  24  through  28 
for  •>  -  9  dB.  An  inspection  of  these  curves  shows  that  the  S/N  required  for  detection 
increases  rapidly  as  •>  increases.  In  fact  the  S/N  in  dB  appears  to  increase  linearly  as 
■i  increases,  doubling  as  ..  goes  from  3  to  6  dB  and  increasing  by  a  factor  of  3  as  <>  goes 
from  3  to  9  dB. 

It  is  interesting  to  compare  the  set  of  curves  for  >  =  3  dB  and  N  1  (no  integration) 
with  those  obtained  using  a  Gaussian  model.  Skolnik  (5)  shows  that  the  S/N  required  com¬ 
pares  favorably  with  the  results  using  a  log-normal  model  (Fig.  14)  as  should  be  expected 
lor  the  >i  =  3  dB  case. 
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normalized  envelope  amplitude  of  signal  plus  cl-  .ter,  l  =  v/v„, 


Fig.  7  -  Probability  densities  of  a  signal  plus  log-normal 
clutter  for  a  -  3  dB,  voltage  S/N  =  2,  and  number  of  pulses 
N  ■=  1  l,  and  3.  Note  aa  N  increases  the  curve  loo, lS  more 
like  a  normal  distribution. 


Fig.  8-  Bias  levels  for  false  alarm  probabilities 
from  10-2  to  1U'H  for  N  -  3  with  <>  -  3,  6,  and 
9  dB 
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Fig.  13  -  Bias  levels  for  false  alarm  proba¬ 
bilities  with  the  number  of  pulses  integrated, 
N  =  1,  3,  10,  and  30  for  a  -  9  dB 


Fig.  14  -  Probability  of  detection  vs  S/N 
for  PfA  =  10'2  to  10’8  with  a  =  3  dB  and 
N  -  1 
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Fl«.  .iO  -  Probability  of  detection  va  S/N  for 
I'...  l0-i  to  10"H  with  •!  a  (i  dB  and  N  =  3 
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Fig.  <21  -  Probability  of  detection  vs  S/N 
for  rFA  =  10-2  to  10‘8  with  n  =  6  dB  and 
/v  =  10 
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Fig.  38  -  Probability  of  detection  vs  S/N 
for  N  -  1,3,  10,  and  30  with  «  =  9  dB  and 
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